Covariant boost and structure functions of baryons in Gross-Neveu models 
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Baryons in the large N limit of two-dimensional Gross-Neveu models are reconsidered. The time- 
dependent Dirac-Hartree-Fock approach is used to boost a baryon to any inertial frame and shown 
to yield the covariant energy-momentum relation. Momentum distributions are computed exactly 
in arbitrary frames and used to interpolate between the rest frame and the infinite momentum 
frame, where they are related to structure functions. Effects from the Dirac sea depend sensitively 
on the occupation fraction of the valence level and the bare fermion mass and do not vanish at 
infinite momentum. In the case of the kink baryon, they even lead to divergent quark and antiquark 
structure functions at x = 0. 
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I. INTRODUCTION 

Given the notorious difficulty of solving non- 
pcrturbative problems in quantum field theories (QFT) 
systematically, it is of theoretical interest to study ex- 
actly soluble model field theories. From the point of 
view of strong interaction physics, a particularly grati- 
fying example is the Gross-Neveu model family in 1+1 
dimensions with Lagrangian [l[ 

C = i> - m ) </> + y [ W) 2 + A(^i 75 ^) 2 ] ■ (1) 

It describes N species of self-interacting fermions (fla- 
vor labels are suppressed as usual) with discrete (A = 0, 
Gross-Neveu (GN) model) or continuous (A = 1, Nambu 
Jona-Lasinio (NJL) model Q) chiral symmetry, possi- 
bly broken by the bare mass term ~ mo- Throughout 
this work we are only interested in the 't Hooft limit 
N — > oo, Ng 2 = const. 

The simple Lagrangian ([1} gives rise to a number of in- 
teresting phenomena such as asymptotic freedom, spon- 
taneous symmetry breaking, hadronic bound states and 
non-trivial phase diagrams as a function of temperature 
and chemical potential. These non-perturbative phenom- 
ena are either accessible through analytical methods, or 
else numerically to any desired accuracy. So far, the main 
focus has been on static properties such as the hadron 
spectrum pj-fllT] or equilibrium thermodynamics ■ 
We see no reason why such models should not be in- 
structive for dynamical problems as well. In principle, 
it is known how to generalize the semiclassical methods 
adequate in the large TV limit to the time-dependent case 
[1, [Tt], [l8|]. For baryons in particular, this amounts to 
replace the relativistic Hartree-Fock (HF) approach by 
its time-dependent generalization. The non-relativistic 
version of time dependent Hartree-Fock (TDHF), intro- 
duced originally by Dirac fl9| . is well-known in nuclear 
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and heavy- ion physics [20|. This paper is a first step to- 
wards attacking time dependent problems in QFT with 
relativistic many-body techniques. We hope to extend 
thereby the spectrum of questions which can be solved 
exactly in field theory models, serving as a testing ground 
for other approaches. 

Consider a baryon in the Gross-Neveu model, i.e., a 
HF state of n valence quarks with the polarized Dirac 
sea. Dashen, Hasslacher and Neveu (DHN) have already 
addressed the issue of covariance of such bound states 
formally from the point of view of semiclassical methods 
and concluded that the relativistic energy-momentum re- 
lation E = V M 2 + P 2 is satisfied [l7|. Here we take up 
this question using relativistic many-body methods. As- 
suming that the baryon state has been constructed in 
the rest frame, we boost it to an arbitrary inertial frame 
and obtain a self-consistent solution of the Dirac-TDHF 
equations. We then compute energy and momentum of 
the moving baryon and verify covariance in the large N 
limit. In general, to compute the energy and the mo- 
mentum of such bound states in QFT is rather delicate. 
It requires renormalization, vacuum subtraction and the 
use of a non-covariant cutoff and a finite box at interme- 
diate stages, at least in the previous approaches in the 
rest frame E3|- In view of these complications, we 
find it worthwhile to examine in detail how covariance 
is maintained, even if the outcome is known on general 
grounds. As a matter of fact, we will not need the ex- 
plicit solution of the baryon in the rest frame (which is 
known analytically only for the GN model, but not for 
the NJL model). This sheds some light on how covari- 
ance emerges in the canonical formalism where it is often 
deeply hidden. 

Having confirmed covariance, as an application we 
compute quark momentum distributions for the baryon 
moving with arbitrary velocity. Here we do need the ex- 
plicit HF wave functions, so that we can carry out an 
analytical calculation for the DHN baryons of the GN 
model only. By taking the infinite momentum frame 
limit, we also determine structure functions, the quan- 
tities of central interest in deep inelastic scattering and 
quantum chromodynamics (QCD), see e.g. [2l|. In such 



2 



a model study it is easy to disentangle valence from sea 
quarks and to study explicitly the dependence of the 
structure functions on the occupation of the valence level 
(or baryon number) and quark masses. 

We will also address the question to which extent 
the HF calculation simplifies in the infinite momentum 
frame, relevant for attempts to work in light-cone quan- 
tization [22L I23I ] . In the meson problem (quark-antiquark 
bound state), such simplifications do indeed occur and 
have been instrumental for the first successful applica- 
tion of light-cone quantization by 't Hooft to large N 
QCD2 [24j. Here as well as in the GN model, equiv- 
alence of the meson spectrum with equal time quanti- 
zation has been firmly established long ago |25l - l28l |. In 
the case of baryons, the situation is less clear. We are 
aware of a few early attempts to compute baryon masses 
on the light-cone for finite N QCD2 [III [3(| and of the 
works 31-l33j relevant for the large N limit, but there 
has been no evaluation of the DHN baryons on the light- 
cone to the best of our knowledge. In the meson case, a 
pure valence quark approach becomes exact in light-cone 
quantization, technically reducing the relativistic random 
phase approximation to the simpler Tamm-Dancoff ap- 
proximation. By contrast, the role of the Dirac sea for 
baryons in the infinite momentum frame is still unclear. 

This paper is organized as follows: In SecQll we briefly 
recall the evaluation of the baryon mass in the rest frame 
and show how this calculation can be simplified signifi- 
cantly. This explains several cancellations which were 
observed in earlier works and yields a new and surpris- 
ingly simple relationship between the baryon mass and 
the self-consistent HF potential. This step is necessary to 
keep the proof of covariance in later sections tractable. In 
Sect. Hill we set up the Dirac-TDHF approach for boost- 
ing the baryon. Sect. [TV] is dedicated to the calculation 
of baryon energy and momentum in an arbitrary inertial 
frame and to the test of covariance. In Sect. [V] we define 
and compute quark and antiquark momentum distribu- 
tions for the moving DHN baryon of the GN model with 
(broken or unbroken) discrete chiral symmetry. Upon 
performing the limit to the infinite momentum frame, 
we arrive at simple expressions for structure functions 
in Sect. IVIl We also discuss issues like valence and sea 
quark contributions, or the dependence on the baryon 
number and the bare fermion mass. We finish this paper 
with a concluding section, Sec. IVHI 



II. SIMPLIFYING THE CALCULATION OF 
THE BARYON MASS IN THE REST FRAME 

The calculation of the baryon mass in the GN model 
Jl @ or of the Shei bound state in the NJL model H, 
|l0( is rather involved. The corresponding calculation of 
energy and momentum of the moving baryon is even more 
difficult. If one looks at the HF calculation in the rest 
frame in detail, one finds that many terms cancel at the 
end. This points to the presence of simplifying features. 



In order to take advantage of these simplifications in the 
case of the moving baryon, we first have to gain a better 
understanding of the calculation in the rest frame. This 
is the topic of the present section. 

Assume that we have solved the HF problem self- 
consistently for given baryon number, 

[-ijrA + j°S(x) + i7 X AP(a:)] <f, a (x) = E a <t> a (x), (2) 

with scalar and pseudo-scalar potentials 

occ 

S(x) = -g 2 y^^4> a (x)4> a (x) +m , 

Q 
OCC 

P(x) = -g 2 ^2(j) a (x)i^5(t)a(x). (3) 

a 

The sums run over all occupied states, i.e., the valence 
level and the Dirac sea. The Hamiltonian and momentum 
operator of the GN model read 



H 



+A(^i 75 ^) 2 ] } , 
P = [dx (-i^d x i>) . 



(4) 



Let us write the expectation values in the HF baryon 
state as 



(H) = J dx(H(x)), 
(P) = J dx(V(x)). (5) 
where the energy and momentum densities are given by 

OCC s ^ 

(H) = \~ i^Te^a + ^{m + S)(f> a (j> a 



(0) 



+-AP0 Q i750 

a ( 1 



We have to subtract the vacuum contribution from {%). 
For ease of notation, we denote the expectation values 
(H),(P} generically by J- and their contribution from 
single particle state a by F a . 



(7) 



The c- number densities T a are x dependent, approaching 
a constant T%. y for x — > ±00. After splitting off the 
constant term, 

= {J~ a — -^asy) + -^asy := -^Sc + -^asy; (8) 
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the localized part J~\ oc can be integrated over x from — oo 
to +00 with a finite result. The constant, asymptotic 
part .Fasy requires the usual careful vacuum subtraction. 
To this end we enclose the system temporarily in a box of 
length L and impose periodic boundary conditions. F^ sy 
is non- vanishing only for the continuum states and differs 
from the corresponding vacuum quantity due to a change 
in the density of states, 



■pa 

asy 



•F asy(^) — -F vac(fc) I 1 



1 dS(k)\ 
L dk ) 



(9) 



Here, k labels the asymptotic momentum of the scatter- 
ing state a and S(k) is the scattering phase shift. Owing 
to the finite box, the fcrmion momenta arc discrctized as 
fe„ = 2irn/L in the vacuum and as 



with the Lagrangian density C from Eq. ([T]). Working 
out its components, we find 



1111 



= H 

g 2 - 



Tn = U - moV'V': 
T01 = i^d x il> = -V, 

T10 = -V-'-d^, 



(15) 



where the divergence of the axial current jg = ip^^^tp 
is given by 



0^ = 2(m - (1 - A) 5 2 ^)^75^. 



(16) 



k n ^S(k n ) 



(10) 



in the baryon state. Subtracting the vacuum contribu- 
tion, we evaluate 



(•^asy(^n) -7"vac(&°)) 



E 



3~&sy{k n ) 



n \ 




2tt k 



lim [5(k)T vac (k) - 6(-k)F vac (-k)} 



(11) 



where we have used Eq. © and taken the limit L — > oo. 
Hence the vacuum subtracted asymptotic part is a pure 
surface term in momentum space. The vacuum energy 
and momentum densities are 



£(^vac(fc)) = - 

L(V vac (k)) = kN 



2fc 2 



mrriQ 



2E k 



N. 



(12) 



Anticipating that S(k) ~ 1/fc for large |fc| (as will be con- 
firmed in the cikonal approximation below), the asymp- 
totic, cc-independent contribution to the baryon energy 
and momentum in the rest frame finally becomes 



Lj2[(^y(k n ))-(n, 



: (A£))] = - lim kS(k), 

7T k— >oo 



0. 



(13) 



We now turn to the localized part. Here it is useful to go 
back to the energy momentum tensor, since its local con- 
servation law leads to drastic simplifications in 1+1 di- 
mensions. The canonical energy momentum tensor reads 

M 



For expectation values in any stationary state, the con- 
servation law 



reduces to 



d»%„ = 

d x {% - m tptp) = 0, 
= 0. 



d x ( V + -d xP 



(17) 

(18) 
(19) 



Here, we have made use of the elementary fact that (j^ 



7° - 7 1 

J5 — J ' 



3l=f=P 



(20) 



in 1+1 dimensions. For the local pieces, we therefore get 

(Hloc) = mo^VOioc, 

(Pioc) = ~\d x p. (21) 

The first equation implies in particular that the energy 
density in the chiral limit of the GN model should be- 
come x-independent. This is responsible for cancella- 
tions encountered if one evaluates all terms separately, 
see the GN model and the Shei bound state 0, [re- 
integrating the densities ([21]) over x and adding the x- 
independent, vacuum subtracted contributions (|13[) . we 
arrive at a simple expression for the baryon mass 



N 



lim kS(k)+m / dx(ipip)i oc , (22) 

7T k— >oo 



and trivially Pb = in the rest frame. It is now easy to 
express the right-hand-side of Eq. (|2"2"]l in terms of the HF 
potential. The asymptotic behavior of the phase shift can 
be determined with the help of the eikonal approxima- 
tion (Glauber theory [IH). To this end, we write down 
the stationary Dirac equation, eliminate the lower spinor 
component v and use the following ansatz for the upper 
component, 



(14) 



u(x) = u(x)e 



ikx 



(23) 
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with slowly varying modulation u[x). In the high energy 
limit one then finds a first order differential equation 

^ = 4 (S 1 + iXP' -S 2 - XP 2 + m 2 ) u (24) 
ax 2k 

from which the asymptotic scattering phase shift for 
Ifcl — > oo can be deduced as 



S(k) 



1 

2k 



/OO 
dx (S 2 + XP 2 
-oo 



(25) 



(assuming no surface term from the derivatives). The 
chiral condensate on the other hand is directly related to 
the HF potential by self-consistency, 

f°° - mjN f°° 
m / dx(4>-ip)i oc = / dx(S - m). (26) 



(27) 



We have introduced the confinement parameter [3j: 



7 



7r mo 
Ng 2 m 



The final relation between the baryon mass and the self 
consistent potential is then 



M B 



N f°° 



+ 7?n) 2 + AP 2 - (m + 7?7i) : 



(28) 

This formula holds in massless and massive GN or NJL 
models. If we insert the known potential S for the mas- 
sive GN model for instance, we recover the result given 
in Eq. (|86|) below. We have also checked the formula 
for the baryon of the massive NJL model in the deriva- 
tive expansion Q and found agreement with known re- 
sults up to 0(m^ 1 ). In the chiral limit of the NJL model 
(A = 1, 7 — 0\ one sees nicely the appearance of a mass- 
less baryon [37| if the potential traces out the chiral circle, 
S 2 + P 2 = m 2 . 

Concluding this section, we note that the cancella- 
tions observed in previous calculations of the baryon 
mass have two distinct sources: The fact that the vac- 
uum subtracted, constant part of the energy density is a 
pure surface term in momentum space, and local energy- 
momentum conservation relating the x-dependent part 
of the energy density to the subtracted chiral conden- 
sate. These observations will help us to evaluate energy 
and momentum of the moving baryon more efficiently in 
Sect. HVl 



III. 



TIME DEPENDENT HARTREE-FOCK AND 
THE BOOSTED BARYON 



Let S(x),P(x) denote the self-consistent scalar and 
pscudoscalar potentials for the baryon in the rest frame. 
The Dirac-TDHF equation in a frame where the baryon 
is moving with velocity v reads 

[vfd^ - S(j(x - vt)) - i7 5 AP(7(x - vt))] ip a (x,t) = 

(29) 



with 7 = (1 — v 2 ) x / 2 (not to be confused with the con- 
finement parameter) and the self-consistency conditions 



S(-y(x — vt)) 



g 2 ^ip a (x,t)ip a (x,t) + m , 



P(~/(x-vt)) = -g 2 Yy a (x,t)i-/Mx,t). (30) 



It is straightforward to solve Eqs. (|29I30[) by a Lorentz 
boost from the rest frame of the baryon to the frame in 
which it has velocity v. Starting point is the ansatz 



il> a {x,t) =Af a e^ 5 / 2 e- iE °< t '(t) a (x') 



(31) 



where £ = artanhw is the rapidity parametrizing the 
boost 



cosh£ — sinh£ 
-sinh£ cosh£ 



(32) 



and E a , <j) a denote the eigenvalues and eigenspinors of 
the HF Hamiltonian in the rest frame. Af a is a (real) 
normalization factor. Insert Eq. (|3"Tj) into Eq. (f2l)|) . 

[i 7 °d t + i 7 1 d x ~ S(x') ~ i 75 AiV)] c^/V^'V^z') = 

(33) 

Next we pull the spinor boost matrix exp(^7s/2) through 
the 7 M matrices to the left and divide it out, using 



-75C/2 ( 7° ^ o7 5(72 _ ( cosh£ sinh£ 
sinh £ cosh £ 



and 

di 
d x 

This yields 



cosh£ — sinh£ 
- sinh £ cosh £ 



Or 
d x < 



(34) 



(35) 



[rffy + i7^ - S(x') - i-f 5 XP(x')] e"^ 4 <j> a {x') = 

. (36) 

or, equivalcntly, the stationary Dirac-HF equation in the 
rest frame, 

[-^d x +fS{x) +i 1 1 XP(x)] 4> a {x) = E a 4> a {x). (37) 

The self-consistency in the boosted frame follows from 
the assumed self-consistency in the rest frame since 

ip a (x,t)ip a (x,t) = A/24> a (x')(f> a (x'), 
ii a {x,t)vy 5 ijj a {x,t) = 7V 2 a (a;')i75'/ ) Q ! (a;')- ( 38 ) 

The normalization factor M a is needed so that the sums 
over occupied states in Eqs. (|30[) arc done correctly in 
both frames of reference. 

In summary, once one has solved the HF equations 
in the rest frame of the baryon, a standard kinematical 
boost is sufficient to transform the solution into a solution 
of the TDHF equations (|29I30|) . All we have to do in the 



following is to compute observables using the boosted 
wave functions, Eq. (f3"Tj) . 
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IV. TEST OF COVARIANCE VIA 
CALCULATION OF ENERGY AND 
MOMENTUM OF THE BOOSTED BARYON 

Let us evaluate the expectation value of the operators 
H and P, Eqs. (j4|, for a moving baryon. Eqs. (|5I6|) 
remain valid provided we replace the rest frame single 
particle spinors <$> a (x) by the boosted ones, ip a (x',t'), 
and S(x),P(x) by the boosted potentials S(x'), P(x'). 
We consider first the contributions (Ha), (P a ) from one 
single particle state a to the sum in Eqs. (|6]). Using 



d x = cosh^Sj;' — sinli£<9t< , 
c ?7s _ cosh £ + 75 sinh £, 



(39) 



and eliminating the time derivative of tp a with the help 
of the Dirac-HF equation 

id t ,ip a (x',t') = [-i l5 d x/ + 1 a S(x') + i 1 1 XP(x')] i> a (x',t') 

(40) 

a large number of terms is generated which we organize 
(in anticipation of the result) as follows, 



1=1 

CPa) = A/2£>«>, (41) 



i = l 



with 



= cosh(20^(a; , )75-^</) a ( a; ') 
i 

+i [m + cosh(2£)S(x')] ^ a (x')K(x') 
+- cosh(2^)AP(x')^a(a; / )i750 

= sinh(2C)4(.T')-^0 a (x'), 
l 

(^i 3) ) = ^sinh(20 [5(x')0a(x')i750a(x') 

-AP(z')<Mz')<M*')] , 

(P«) = smh(20<i>l(x'h 5 -d x >Mx') 

l 

+ isinh(2O5(x')^(x')0a(x') 

+ - S mh(2$)\P(x')0a(x')h5(l>a(x'), 

(V&) = co S h(204>i(x')-d x ,Mx'), 

1 

(V a 3) ) = isinh 2 e[5( 2; , )^(a ; , )i750a(x') 

-AP(a/)^ a (af')^a(s / )] • ( 42 ) 

Once again we have to treat separately the localized, sub- 
tracted parts of the various densities and their constant, 
asymptotic parts. Consider the localized densities first. 
The normalization factor Af a is necessary to transform 



the sum over continuum states from one frame into the 
other one. To leading order in 1/L needed here, it is 
given by 



N 2 = 

J v rv 



E 



(43) 



We denote the center-of-mass (CM) frame kinematical 
variables by (E, k), the laboratory (LAB) frame variables 
by (w, q), so that 



cosh £ sinh £ 
sinh £ cosh £ 



(44) 



When integrating over q, Lorentz invariance of the mea- 
sure dq/u allows us to relate the quantities summed over 
all continuum states in the LAB and CM frames. For the 
discrete states, N a = 1. The 2nd issue is the integration 
over x when going from densities to expectation values. 
If we work in a finite box of length L and transform in- 
tegration variables from x to x', the integration limits 
acquire a time dependence, 

L/2 rL/2 cosh £—t sinh £ 

da; -> — / dx'. (45) 

-L/2 c °sh £ J-L/2 cosh £-t sinh £ 

For the localized densities, we may safely extend the inte- 
gration limits to ±oo since we are eventually interested in 
the limit L — > oo (for finite L, we assume that the times 
considered are such that the baryon does not yet see the 
walls). It is then clear that the i = 2 pieces in Eq. (|4"Tj) 
do not contribute to the expectation values of H, P, be- 
ing proportional to the baryon momentum in the rest 
frame. Moreover, the i = 3 pieces vanish due to parity 
(the integrand is odd under reflection). The only non- 
vanishing contributions left are the i = 1 terms. They 
simplify drastically even before integrating over x' once 
we invoke the local conservation of the energy momen- 
tum tensor, Eq. (|21[) . As a result, the contribution from 
the localized densities to the expectation values of H, P 
is simply 



(H) ioc = cosh£m / dx(ipip)i oc , 

J — OO 
p OO 

(P)ioc = sinh£m / dx(-0^)i oc . (46) 



Hence this part is covariant by itself, cf. the correspond- 
ing contribution to the baryon mass in Eq. (|22[) . 

Now consider the constant, asymptotic terms in the 
densities. The expectation values for the continuum state 
labelled by k in the CM frame and by q in the LAB frame 
are given by [cf. Eqs. (|9I12I) ] 

AT 2g 2 + m 2 +mm Q ( 1 d6(k(q)) \ 
L(U a ) &SY N 

/ l d6(k(q)) \ 
asy = Nq\l — — — — — ] . (47) 



L(V a , 



L dq J 
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If F denotes temporarily H or P and T the corresponding 
densities W, V, the relation analoguous to Eq. (ITT1) for the 
boosted baryon reads 



(F) 



asv 



L f A/2 P-^^Him^M)} 

J -A/2 27rd< 7 



= -— lim [S(k(q))(T mc (q)} 
-5(k(-q))(T vac (-q))]. 



(48) 



As shown above in Eq. (|25p . for large k the phase shift 
behaves as 



*< 4 H +o (j!» 



where 



i] = lim = — / dx (S 2 + XP 

fe^oo 2 /.^ 



2 - m 2 ) 



and hence 



q 



for q — > ±oo 



(49) 



(50) 



(51) 



for negative energy states (use k = cosh£g — sinh£w) 
Together with the trivial free asymptotic behavior 



L(H vac (q)) « -\q\N 
L(V vac (q)) « giV 



(52) 



for q — > ±oo, this yields the following final result for the 
contribution from the asymptotic densities to energy and 
momentum of the boosted baryon, 



N 

(H)asy = cosh£ — lim kS(k), 

7T k— >oo 

N 

(-P)asy = sinh£ — lim kS(k). 

7T k— >oo 



(53) 



Eqs. (|22l46l53p then confirm the covariance of the spec- 
trum, 



(H) = cosh ^ M B , 
(P) = sinh£M B , 



(54) 



even without invoking the explicit solution of the GN 
baryon. This is important since such a solution is known 
analytically only for the GN model (A = 0), but not for 
the massive NJL model (A = 1). 



V. QUARK AND ANTIQUARK DISTRIBUTION 
FUNCTIONS FOR THE BOOSTED BARYON 



In the preceding sections, we have verified that the 
TDHF equation for a baryon moving with velocity v 
can be solved by a Lorentz boost of the HF solution 
in the rest frame. The covariant energy-momentum re- 
lation was found in such an approach. This supports 



strongly that the Dirac-TDHF approach is the correct 
procedure in the large N limit. In general, being able 
to work in different Lorentz frames does not yet yield 
any new physics insights, even if it is satisfying from a 
theoretical point of view. Nevertheless, there is one in- 
teresting application which we shall pursue in this and 
the following section: The study of distribution func- 
tions for quarks and antiquarks in momentum space as a 
function of baryon velocity. By going to the infinite mo- 
mentum frame limit, we can then extract structure func- 
tions closely related to those appearing in QCD analyses 
of dccp-inclatic electron-proton scattering. Previous at- 
tempts to determine structure functions in various quark 
models (e.g., the MIT bag model) have suffered from a 
lack of covariance of the underlying formalism, a prob- 
lem which we do not share. On the other hand, structure 
functions are exceedingly hard to compute in lattice QCD 
where one is typically able to evaluate the lowest few mo- 
ments only (see e.g. HHHil). Therefore it may be useful 
to have a covariant field theoretic model where one can 
compute structure functions exactly. For this purpose, 
we need the explicit HF wave functions of the baryon 
which arc known analytically only for the massless and 
massive GN models [A = in Eq. (TTJ)]. In order to get as 
much analytical insight as possible, we restrict ourselves 
to these variants of the GN model family from here on 
and deal exclusively with DHN baryons @, H, § . 

The fermion field operator can be expanded either in 



terms of the free basis, 



or the HF basis, ipn 



(±) 



k 

= £ (A n 4+^ (x) + B n 4-\x)) . (55) 



The superscript (±) refers to positive and negative energy 
states. The annihilation operators in the two bases are 
related by the Bogoliubov transformation 

n 

h = ^{(^Ww^ + ^^i-'R}. (56) 

n 

We define (momentum space) distribution functions for 
positive and negative energy fermions as 



(HF|4a fe |HF) = £|(4 + \4 +) )| 2 +El(^ +) >^ 



2 

n J I 3 



(HF|&i& fe |HF) = £|(^-\^ +) )| 2 ±]T|(</4-W 



*H 4-))| 2 , 



(57) 

where |HF) denotes the HF baryon state and the sums 
run over all occupied states, i.e., the filled Dirac sea and 
the valence level. In view of the physical interpretation 
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we convert the negative energy fermions into antiparticlcs 
by means of a standard particle-hole conjugation, 



b\b k = d^ k d f _ k = 1 - d/_ k d-k, 



(58) 



and identify momentum distribution functions for quarks 
and antiquarks as follows, 



W q (k) = (HF|4o fc |HF), 
W q (k) = l-(HF|6i fc 6_ fc |HF). 



(59) 



Note that there is no summation over flavor, so that all 
momentum distributions refer to a single flavor. 

For DHN baryons and antibaryons (fcrmion number 
±vN), we display the discrete and continuum states ex- 
plicitly, 



W q B (k) = v\(<t>l + \4 +) )\ 2 + \(<t>l +, ^o~>)\ 

cont 

(+) ./.(-hi* 



(+) ,d-)\fl 



wm) = i-H(^ fc \4 +) )iM(&v>ni 



/,(-) »/,(-)m2 



cont 



cont 

wf(fe) = l-a-^K^,^)! 2 -^!^'^)! 2 . 

(60) 

Remember that the GN model is charge conjugation sym- 
metric, so that all single particle states come in pairs 
with opposite energy. Eqs. ([60j) merely reflect the dif- 
ferent ways in which the discrete levels are filled. In the 
baryon, the negative energy state is completely filled with 
TV fermions, the positive energy state partially with fill- 
ing fraction v = n/N. In the antibaryon, the positive 
energy state is empty whereas the negative energy state 
is filled with fraction 1 — v. By charge conjugation (C), 
we must have 



W q B (k) = W B (k), W q B (k) = W B {k). 



(61) 



This enables us to express the C-odd combination in 
terms of valence quantities only and to simplify some- 
what the C-even combination, 

W q B (k)-W B (k) = v{\(<f>£\l>P ) )\* + \tt?\4~ ) )\*}> 

w B (k) + w B (k) = H(4 +) ,4 +) )l 2 

+(2-i/)|(4 +) ,^" ) )| a 



(62) 



Since the single particle wave functions for the DHN 
baryon are known, the distribution functions for the dis- 
crete states can easily be evaluated analytically. We refer 
the reader to @ for the relevant detailed wave functions 
in the rest frame of the baryon. Upon boosting these 
spinors according to Eq. (|31[) and using units such that 
the dynamical fermion mass in the vacuum has the value 
m = 1 from now on, we find 



L 
2^ 



i(*i + \4 +) )i a 



with 



7T (cos(aA_) — qsin(aA_) + E q ) 
8E k y cosh 2 £ cosh 2 (/3AJ) 
7r (cos(oA + ) + gsin(aA + ) — E q ) 



8E k y cosh 2 £ cosh 2 (/3A+) 



(63) 



2c 



y cosh £ ' 

7T 

2y cosh£ ' 



Co 



-artanh y 



A± = y/l - y 2 sinh£± k, 
and the boosted kincmatical variables 



E„ 



i - 



cosh£ sinh£ \ / E k 
sinh£ cosh£ J \ k 



(64) 



(65) 



The factors L/(2tt) in Eq. (62) have been introduced with 
regard to the limit L — >• oo. The parameter y is deter- 
mined by the occupation fraction v of the positive energy 
valence state of the DHN baryon [3j , 



/TTV 

V = sin 



(66) 



Note also the symmetry relation which can be used to 
obtain other related matrix elements, 

\(4 a \4 v) )\ 2 = \(^7\4' v) )\ 2 (<™ = ±i). (67) 

Matrix elements involving continuum HF states have 
been computed as follows: We start from box normal- 
ized spinors and boost the HF spinors [f| from the rest 
frame to the moving frame. The relevant matrix element 



Scant =(£) U^^K)? 



(68) 



with the analytical result 

_ (K - g) siii(o-C) - (gK + 1) cos(aC) + E q E K 



8 cosh 2 £(K 2 + y 2 )E k E Q sinh 2 (/3C) 



Here. 



Eq 



cosh £ — sinh S, \ f Ek 
— sinh £ cosh £ I \ K 



(69) 
(70) 
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a and /3 are defined as in Eq. (|64[) . and 
C = k-Q. 



(71) 



In the course of this calculation, one runs into the Fourier 
transform of tanhz with finite integration limits. This 
was done as follows: Add and subtract the function sgnz. 
In the Fourier transform of (tanh z— sgn z), we may safely 
extend the integration limits to ±00 in the large L limit, 



dz(tanhz — sgnz)e~ 



l))z 



-2i 



tt/2 



1 



sinh(7rp/2) p 



. (?2) 

The integral over sgnz is then evaluated with finite in- 
tegration limits. Performing the limit L — s- 00 at the 
end, all momenta go over into continuous variables and 
we normalize the momentum distributions accordingly. 
Summarizing, our final result for the discrete and contin- 
uum contributions is given in closed analytical form up 
to a one-dimensional integration, 



W q B (k) 
W B {k) 

val(^) 



W B 



vWi(k) + W 2 (k) + W 3 (k), 
(l-u)W 2 {k) + W 3 (k), 
vW x {k), 



(73) 



with 



Wi(k) = 

W 2 (k) = 
W 3 (k) = 
dQ 



tt (cos(aA_) — qsin(aA_) + E q ) 
8E k y cosh 2 £ cosh 2 (/3A_) ' 
7r (cos(aA + ) + gsin(aA + ) — E q ) 
8E k y cosh 2 £ cosh 2 (ft A+) 



(K - q) sin(aC) - (qK + 1) cos(aC) + E q E K 



! cosh 2 £(K 2 + y 2 )E k E Q sinh 2 (/3C) 



(74) 



These distribution functions are normalized according to 
the following "sum rules" to baryon number and baryon 
momentum, 



n 
N 

Pb 

N 



dk (W q B (k)-W q B (k)), 
dkk(W B (k) + W B (k)) . (75) 



Let us illustrate the behavior of the distribution func- 
tions with a few examples for DHN baryons in the mass- 
less GN model. We first choose a moderate occupation 
fraction of the valence level {v = 0.75). In Fig.Q] we show 
the quark and antiquark distribution functions W B and 
W B for a baryon at rest, together with the contribution 
from the positive energy discrete state, W B al . Fig.[5]is the 
corresponding plot for baryon momentum Pb/N = 5 (in 
units where m = 1). At this value of v, the contributions 
from the Dirac sea and from antiquarks are rather small, 
with little dependence on baryon momentum. This pic- 
ture changes if we go into the regime where the bound 



15- 



/0.4- 



0.3 



0.2 



0.1- 



-2 



k 



FIG. 1: Quark (W q , thick curve), antiquark (Wq , thin curve) 
and valence quark (W^, dotted curve) distribution functions 
for baryon at rest versus fermion momentum. The parameters 
are v = 0.75, 7 = 0, units m = 1. 




FIG. 2: Same as Fig. [TJ but for a boosted baryon with mo- 
mentum Pb/N — 5. Note the changed scale on both axes. 



state is highly relativistic by choosing v = 0.99999. As 
shown in Figs. |3]and|4l here the valence quark distribu- 
tion function bears little resemblance with the full quark 
result, showing an enhanced role of the Dirac sea. To 
quantify the dependence of sea effects on total baryon 
momentum, we have also integrated the quark and anti- 
quark distribution functions over all momenta k. At this 
point it is convenient to introduce hatted quantities for 
baryon momentum, baryon mass and fermion number 



N ' 



M B = 



Mb 
N ' 



N 



(76) 



to get rid of trivial N- dependences. In Figs. [5] and [6] 
are shown the contributions from quarks and antiquarks 
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FIG. 3: Same as Fig.Q] baryon at rest but different occupation 
fraction. Parameters v — 0.99999, 7 = 0. 
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FIG. 4: Same as Fig. [3] but for boosted baryon with momen- 
tum Pb/N = 5 and adjusted scale. 



to reduced fermion number versus reduced baryon mo- 
mentum. At v = 0.75 (Fig. [5]), sea effects are smafl ev- 
erywhere, quickly reaching some (non-zero) asymptotic 
value. At v = 0.99999 (Fig. U), they drop rapidly be- 
tween Pb — and Pg = 10, but then stay constant at 
a sizeable level. The lesson we draw from this is that 
antiquark effects apparently do not disappear in the in- 
finite momentum frame, although they are reduced sig- 
nificantly as compared to the baryon at rest. This is in 
contrast to the earlier observation that antiparticle ef- 
fects in the structure of mesons are completely quenched 
in the infinite momentum frame [40]. It suggests that a 
light-cone approach to baryons may be less efficient than 
for mesons, at least in the large N limit. 



2- 



Nt 



-4- 



10 



Pb 



20 



25 



30 



FIG. 5: Dashed straight line: reduced fermion number Nf = 
v = 0.75 of DHN baryon. The solid curves show how fermion 
number is made up from quarks (upper curve, positive contri- 
bution) and antiquarks (lower curve, negative contribution), 
as a function of baryon momentum. The asymptotic values 
reached around Pb = 5 are 0.80 quarks and 0.05 antiquarks, 
as compared to 0.95 quarks and 0.20 antiquarks at Pb = 0. 
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FIG. 6: Same as Fig. but different fermion number Nf = 
v = 0.99999. Here the asymptotic values reached around 
Pb — 10 are 1.55 quarks and 0.55 antiquarks, as compared to 
5.79 quarks and 4.79 antiquarks at Pb = 0. 



VI. INFINITE MOMENTUM FRAME AND 
STRUCTURE FUNCTIONS 



To facilitate the infinite momentum frame limit, we 
scale the momentum variable with the reduced baryon 
momentum, 



xP E 



(77) 
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and introduce the rescaled densities 

wf(x) = P B W q B {xP B ), w B {x)=P B W?{xP B ). (78) 

In this form, the limit £ — > oo can readily be taken. Note 
that since P B is defined without the factor N, x is not 
restricted to [0, 1] in the infinite momentum frame like 
the standard Bjorken variable, but rather to [0, N] (i.e., 
the positive half- axis in the limit N — > oo). The result 
for the positive and negative energy discrete states is 



,(±) 



(x) 



ttM b (1 — sinK^) 
4ycosh 2 (7TK T /(4co)) 



(79) 



with 



k± = 



2c V 1 " V 2 ± m b* 



y 



The negative energy continuum contribution reduces to 



(-) 



(x) 



dqF(q,x) 



(81) 



with 

F(q,x) = 
and 

From Eqs. (|79H83[) . full quark and antiquark structure 
functions and the valence quark structure function for 
the baryon can be obtained as follows, 

w?(x) = vw { d ii(x)+wi;J c {x)+w { c ~l{x), 



M B [{q 2 - 1)cosk - 2gsinK + q 2 + l] 
2 l(q 2 - I) 2 + Vg 2 ] sinh 2 (™/(4c )) 

2c (M B x + q) 

K = . 

y 



(82) 



(83) 



,.(-) 



(*), 



w^ al (x) 



(84) 



They are normalized according to 



1 



dx(wf(x)-w*(x)), 

] 

dxx(wf(x)+w?(x)) 



(85) 



Watch out the integration limits characteristic for the 
large N limit. Finally we recall the relations between 
wand M B holding in the massive Gross-Neveu model 



1,1, 



y 

M B 



- (6 + 7 tan 6) , 

7T 

sin 9, 
2 

— (y + 7 axtanhy) . 



(86) 
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0.5 



1.5 



2.5 



3.5 



(80) FIG. 7: Quark (wf , thick curve), antiquark (wf , thin curve) 
and valence quark (w^, dotted curve) structure functions for 
baryon versus rescaled fermion momentum. The parameters 
are v = 0.75,7 = 0- This graph can be regarded as the 
infinite momentum frame limit of Figs. 11121 with appropriately 
rescaled axes. All structure functions drop to zero at x — 
and vanish identically for x < 0. 




FIG. 8: Same as Fig. [7] for different occupation fraction v — 
0.99999. Continuation of Figs. 13141 to the infinite momentum 
frame. 



In Figs. [7]and[8j we show by way of example structure 
functions for the massless GN model at v — 0.75 and 
v = 0.99999. Apart from trivial rescalings of both axes, 
these plots can be regarded as the continuation of Figs.Q]- 
2] to the infinite momentum frame. We observe that in 
the infinite momentum frame sea effects are important for 
almost complete filling of the valence level. They can be 
reduced by choosing a smaller occupation. An alternative 
way to suppress sea effects is to switch on the bare quark 
mass while keeping the occupation the same, see Fig. |9] 
For the same occupation fraction as in Fig. [5] and the 
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FIG. 9: Same as Fig. [8] but for massive GN model with pa- 
rameters v — 0.99999, 7 — 0.1. Comparison with Fig.[7]shows 
that introducing a bare quark mass has a similar effect as re- 
ducing the occupation fraction, driving the baryon into the 
non-relativistic regime. 




6 0.2 0.4 0^6 o'.& I 

7 



FIG. 10: Dashed curve: reduced fermion number Nf = v = 
1 of baryon in the massive DHN model. The solid curves 
show how fermion number in the infinite momentum frame 
is made up from quarks and antiquarks, as a function of the 
confinement parameter (proportional to the bare quark mass) . 
Both curves show a logarithmic divergence in the chiral limit 
7 — > relevant for the kink baryon, see main text. 
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u 0.1 0.2 0.3 0.4 0.5 

7 

FIG. 11: Decomposition of the total baryon momentum into 
valence quark, sea quark and antiquark contributions, in the 
infinite momentum frame, versus 7. As in Fig. 1101 the valence 
level is fully occupied (y — 1). The two circles on the vertical 
axis (7 = 0) indicate the endpoints of the boundaries separat- 
ing the different contributions and are evaluated analytically 
in the main text. 



the infinite momentum frame. The two curves add up to 
the reduced fermion number 1 and exhibit strong devi- 
ations from a valence picture near 7 = 0. At the point 
v = 1,7 = 0, the parameter y tends to 1 and the DHN 
baryon goes over into a kink-antikink pair at infinite sep- 
aration. Our results indicate that in this limit the baryon 
consists of diverging numbers of quarks and antiquarks, 
even in the infinite momentum frame (below we will ar- 
gue that the divergence is logarithmic). Whereas Fig. [TUl 
details the contributions from quarks and antiquarks to 
the total fermion number, Fig. [Tl] decomposes the total 
baryon momentum in the infinite momentum frame into 
contributions from valence quarks, sea quarks and anti- 
quarks, again at v = 1 as a function of 7. The boundaries 
between the three regions in the plot hit the 7 = axis at 
the points indicated by circles which will be determined 
below. 

The limit y 1 is the ultrarclativistic limit for the 
internal baryon structure. It can be interpreted as refer- 
ring to the kink-like baryon of the massless GN model 
(here with fully occupied valence level). In this limit the 
formulae for the structure functions greatly simplify, 



moderate value 7 = 0.1 of the confinement parameter, 
one gets a picture in qualitative agreement with the one 
at 7 = and v = 0.75 in Fig. 

Finally, we illustrate the influence of the bare quark 
mass on sea effects with the help of integrated quanti- 
ties. Here we consider full occupation [y = 1) and vary 
the confinement parameter 7. Fig. 1101 exhibits the indi- 
vidual contributions from quarks and antiquarks to total 
fermion number as a way of quantifying sea effects in 



lim (x) 



cosh" x Jo 

/>oo 

lim w? (x) = dq 

/n 



dq 



(ir 2 + Aq 2 ) sinh (x + q) 
(87) 



(tt 2 + Aq 2 ) sinh 2 (x + q) 



The 1/cosh 2 term is due to the discrete levels, the in- 
tegral to the negative energy continuum. The sum rules 
for v = 1 arc satisfied. There is an infrared diver- 
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XWIX. 




FIG. 12: Momentum distribution functions xw(x) in the in- 
finite momentum frame and in the kink limit y = 1. Quark 
(xw q (x), thick line), antiquark (xwf(x), thin line) and va- 
lence quark (xw^ a i(x), dotted line) structure functions are 
exhibited. The solid curves hit the x = axes at the point 
2/tv 2 in accordance with Eq. 



gence at x = in the continuum contribution, 

2 

lim w® (x) w lim (x) « — (x — > 0). 



It gives rise to a logarithmic divergence in the number of 
quarks and antiquarks upon integration over x, explain- 
ing the steep rise of the curves in Fig. [101 towards 7 = 0. 
If we plot the momentum distribution xw(x), everything 
is well-behaved and we obtain Fig. which also shows 
the valence level contribution. These are the structure 
functions of the kink baryon. By integrating Eqs. ([57)1 
over x we can evaluate analytically the contribution to 
the total momentum from valence quarks, sea quarks and 
antiquarks. In this way we find 5 In 2 for valence quarks, 
i for sea quarks and | — 5 hi 2 for antiquarks. In other 
words, valence quarks carry 35%, sea quarks 50% and 
antiquarks 15% of the total (kink) baryon momentum in 
the infinite momentum frame. This explains the location 
of the points on the 7 = axes drawn in Fig. 111! 



VII. SUMMARY AND OUTLOOK 

In the past, exactly solvable QFT models have given 
insights mostly into static problems: dynamical masses 
of constituents, the structure and mass of bound states, 
phase diagrams at finite temperature and chemical po- 
tential. All of this touches upon non-perturbative is- 
sues of interest to strong interaction physics. A par- 
ticularly rewarding class of such toy models are large 
N fermionic models in 1+1 dimensions of Gross- Neveu 
type. In the present work, a first attempt was made to 
widen the scope of these studies towards time-dependent 



questions. Due to the fact that lattice gauge theories 
can only be solved in Euclidean time, dynamical non- 
perturbative phenomena are even more elusive in real- 
istic theories than static ones, although they may be of 
considerable interest from the physics point of view. 

At first, we wanted to verify that the large N limit pre- 
serves covariance. Although no problems are expected on 
general grounds, we were interested in a detailed study 
to see whether one can actually boost a composite, rela- 
tivistic bound state in practice. We found that by gener- 
alizing Dirac-HF to its time-dependent version, Dirac- 
TDHF, it is indeed possible to boost baryons to any 
frame and to confirm the covariant relation between en- 
ergy and momentum. This underlines the advantage of 
using QFT toy models rather than more phenomenologi- 
cal quark models, where covariance is always violated at 
some level. We have learned how covariance is restored 
at the end of a calculation involving a non-covariant mo- 
mentum cutoff and a finite box at intermediate stages. 
Renormalizability is a decisive feature here. As a by- 
product, we have understood cancellations observed in 
former calculations of the baryon mass in the rest frame 
and derived a simple formula relating the mass with the 
HF potentials. 

As is well known, internal motion and CM motion do 
not decouple in relativistic bound states, unlike in non- 
relativistic quantum mechanics. In the GN model, one 
can control the extent to which the internal motion is 
relativistic by means of the parameters v (occupation 
fraction of the valence level) and 7 (related to the bare 
mass). Through the choice of velocity v on the other 
hand we can choose the degree to which the CM motion is 
relativistic. Exploring the full range of these parameters 
should enable one to better understand how internal and 
overall motions are interwoven. 

One way of analyzing how the structure of the bound 
state depends on its velocity is through momentum dis- 
tribution of quarks and antiquarks. Since all the wave 
functions are known analytically in the case of the GN 
model with discrete chiral symmetry, such a calculation 
can be done exactly and requires only a one-dimensional 
numerical integration. We were interested in the evo- 
lution of momentum distributions from the rest frame 
(where one has some physical intuition from atoms or 
nuclei) to the infinite momentum frame, which is accessi- 
ble in nature through deep inelastic scattering but where 
one's intuition is in general less developed. We illustrated 
our results with a few examples where relativistic effects 
were either weak or strong, as evidenced by the size of 
sea effects. We had actually hoped to be able to exploit 
covariance in order to simplify the HF problem, similarly 
to what has already been achieved for the meson Bethe- 
Salpetcr equation. Unfortunately, we did not see decou- 
pling of the Dirac sea in the baryon case, although the sea 
effects arc reduced in the infinite momentum frame. This 
sheds light on the possible use of light-cone quantization 
for baryons. 

By performing the limit to the infinite momentum 
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frame, the momentum distributions go over into func- 
tions closely related to structure functions in real QCD. 
They satisfy exact sum rules, but the variable similar 
to Bjorken x cannot be restricted to the interval [0,1] 
but only the positive half axis in the large N limit. Of 
particular interest is the kink limit where the formulae 
simplify further. The ultrarelativistic character of the 
bound state gives rise to a diverging number of quarks 
and antiquarks (coming from the low x region) which has 
little in common with the non-relativistic, valence-level 
type baryon. 

Encouraged by these first results, we plan to ad- 
dress more demanding dynamical problems like kink- 
kink scattering, acceleration of baryons by external fields 
or non-equilibrium thermodynamics in the future. This 



may even be of some interest for the parallel world 
of condensed matter physics where the DHN baryons 
live a life as solitons, polarons and excitons in quasi- 
onedimensional systems, important for example for con- 
duction properties of doped polymers [36 1. 
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